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Abstract 

An approach to the solution of NP-complete problems based on 
quantum computing and chaotic dynamics is proposed. We consider 
the satisfiability problem and argue that the problem, in principle, 
can be solved in polynomial time if we combine the quantum com- 
puter with the chaotic dynamics amplifier based on the logistic map. 
We discuss a possible implementation of such a chaotic quantum com- 
putation by using the atomic quantum computer with quantum gates 
described by the Hartree-Fock equations. In this case, in principle, 
one can build not only standard linear quantum gates but also nonlin- 
ear gates and moreover they obey to Fermi statistics. This new type 
of entaglement related with Fermi statistics can be interesting also for 
quantum communication theory. 
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1 Introduction 



There are important problems such as the napsack problem, the traveling 
salesman problem, the integer programming problem, the subgraph isomor- 
phism problem, the satisfiability problem that have been studied for decades 
and for which all known algorithms have a running time that is exponential 
in the length of the input. These five problems and many other problems be- 
long to the set of NP-complete problems. Any problem that can be solved 
in polynomial time on a nondeterministic Turing machine is polynomially 
transformed to an NP-complete problem [TJ. 

Many NP-complete problems have been identified, and it seems that such 
problems are very difficult and probably exponential. If so, solutions are still 
needed, and in this paper we consider an approach to these problems based 
on quantum computers and chaotic dynamics. 

It is widely believed that quantum computers are more efficient than clas- 
sical computers. In particular Shor [0] gave a remarkable quantum polynomial- 
time algorithm for the factoring problem. However, it is unknown whether 
this problem is NP-complete. 

The computational power of quantum computers has been explored in 
a number of papers. Bernstein and Vasirani || proved that BPPCBQPC 
P SPACE. Here BPP stands for the class of problems efficiently solvable in 
the classical sense, i.e., the class of problems that can be solved in polynomial 
time by probabilistic Turing machines with error probability bounded by 
1/3 for all inputs. The quantum analogue of the class BPP is the class 
BQP which is the class of languages that can be solved in polynomial time 
by quantum Turing machines with error probability bounded by 1/3 for all 
inputs. 

The question whether NPCBQP, i.e., can quantum computers solve NP- 
complete problems in polynomial time, was considered in M. It was proved 
in |4| that relative to an oracle chosen uniformly at random, with probability 
1, the class NP can not be solved on a quantum Turing machine in time 
o(2 n//2 ) . An oracle is a special subroutine call whose invocation only costs 
unit time. This result does not rule out the possibility that NPCBQP but 
it does establish that there is no black-box approach to solving NP-complete 
problems in polynomial time on quantum Turing machines. We would like 
to mention that these results are not immediately applicable to the chaotic 
quantum computer which we consider in this paper. 

For a recent discussion of computational complexity in quantum comput- 
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ing see [^|, [(], [7J. Mathematical features of quantum computing and quantum 
information theory are summarized in [§]. A possibility to exploit nonlin- 
ear quantum mechanics so that the class of problems NP may be solved 
in polynomial time has been considered by Abrams and Lloyd in ||. It 
is mentioned in || that such nonlinearity is purely hypotetical; all known 
experiments confirm the linearity of quantum mechanics. 

The satisfiability problem (SAT), which is NP-complete problem, has 



been considered in quantum computing in |10| . It was shown in |10| that the 
SAT problem can be solved in polynomial time by using a quantum computer 
under the assumption that a special superposition of two orthogonal vectors 
can be physically detected . The problem one has to overcome here is that 
the output of computations could be a very small number and one needs to 
amplify it to a reasonable large quantity 

In this paper we propose that chaotic dynamics plays a constructive role 
in computations. Chaos and quantum decoherence are considered normally 
as the degrading effects which lead to an unwelcome increase of the error 
rate with the input size. However, in this paper we argue that under some 
circumstances chaos can play a constructive role in computer science. In par- 
ticular we propose to combine quantum computer with the chaotic dynamics 
amplifier. We will argue, by using the consideration from [|l(J, that such a 
chaotic quantum computer can solve the SAT problem in polynomial time. 

As a possible specific implementation of chaotic quantum computations 
we discuss the recently proposed atomic quantum computer [TT[. It is pro- 



posed in ||TT| to use a single atom as a quantum computer. One can im- 
plement a single qubit in atom as a one-particle electron state in the self- 
consistent field approximation and multi-qubit states as the corresponding 
multi-electron states represented by the Slater determinant. 

A possible realization of the standard quantum gates in the atomic quan- 



tum computer by using the electron spin resonance has been discussed in |TT 



In this paper we argue that in the atomic quantum computer one can build 
also nonlinear quantum gates because the dynamics of the multi-electron 
atom in the very good approximation is described by nonlinear Hartree-Fock 
equations. 

The tensor product structure of states is very important for computations 
and the multielectron atom admits such a structure. More exactly, instead of 
the standard tensor product used in quantum computing we have to use the 
Slater determinant to take into account the Fermi statistics. The standard 
computational basis in quantum computing does not have Bose or Fermi 



3 



symmetry. In the atomic case we have to make an appropriate modification 
of quantum gates to take into account Fermi statistics and this leads to a 
new type of entanglement related with Fermi statistics. 

Such Fermi or Bose entanglement could be interesting also for quantum 
communication theory, in particular for quantum teleportation [^, [23f1 . 



2 SAT Problem 

Let {xi, ■ ■ ■ ,x n } be a set of Boolean variables, Xi = or 1. Then the set of 
the Boolean variables {xi,x"i, ■ ■ ■ ,x n ,x n } with or without complementation 
is called the set of literals. A formula, which is the product (AND) of dis- 
junctions (OR) of literals is said to be in the product of sums (POS) form. 
For example, the formula 

(xi V x 2 ) (xi) (x 2 V x 3 ) 

is in POS form. The disjunctions {x\ V x 2 ) , (xi) , {x 2 V X3) here are called 
clauses. A formula in POS form is said to be satisfiable if there is an assign- 
ment of values to variables so that the formula has value 1. The preceding 
formula is satisfiable when x± = 0, x 2 = 0, X3 = 0. 

Definition(S AT Problem). The satisfiability problem (SAT) is to deter- 
mine whether or not a formula in POS form is satisfiable. 

The following analytical formulation of SAT problem is useful. We define 
a family of Boolean polynomials f a , indexed by the following data. One a is 
a set 

a = {Si, Sn, Ti, T N ] , 
where Si, Ti C {1, n} , and f a is defined as 

fa{xi, ■ ■ ■ , X n ) = Jl ( 1 + Jl (1 +0 Jl X b J . 

i=l \ aeSi b£T t / 

We assume here the addition modulo 2. The SAT problem now is to 
determine whether or not there exists a value of x = (xi, • • • , x n ) such that 
/«(x) = 1. 
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3 Quantum Algorithm 



We will work in the (n + l)-tuple tensor product Hilbert space 7i = ®" +1 C 2 
with the computational basis 

\x x , ...,x n ,y) = ®" =1 \x t ) ® \y) 

where Xi,...,x n , y = or 1. We denote \x±, x n , y) = |x, y) . The quan- 
tum version of the function /(x) = / Q (x) is given by the unitary operator 
Uf |x, y) = |x, y + /(x)) . We assume that the unitary matrix Uf can be build 



in the polynomial time, see ]I0 |. Now let us use the usual quantum algorithm: 



(i) By using the Fourier transform produce from 1 0,0) the superposition 



(ii) Use the unitary matrix Uf to calculate /(x) : 
\v f ) = U f \v) = ^=J2\^fW) 



2 n 

X 

Now if we measure the last qubit, i.e., apply the projector P = I <g> |1) (1| to 
the state \vj) , then we obtain that the probability to find the result /(x) = 1 
is \\P \vf)\\ = r/2 n where r is the number of roots of the equation /(x) = 1. 
For small r the probability is very small and this means we in fact don't get 
an information about the existence of the solution of the equation f(x) = 1. 
Let us simplify our notations. After the step (ii) the quantum computer will 
be in the state 

\v f ) = yjl-q 2 bo) ® |0> + q \<fi) ® |1> 



where \tpi) and l^o) are normalized n qubit states and q = ^/r/2 n . Effectively 
our problem is reduced to the following 1 qubit problem. We have the state 



-g 2 |0) + g|l) 

and we want to distinguish between the cases q = (i.e. very small q) and 
q > 0. To this end we propose to employ chaotic dynamics. 
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4 Chaotic Dynamics 



Various aspects of classical and quantum chaos have been the subject of 
numerious studies, see [jl2| and ref's therein. The investigation of quantum 
chaos by using quantum computers has been proposed in |L3L [14], [15| . Here 



we will argue that chaos can play a constructive role in computations. 

Chaotic behaviour in a classical system usually is considered as an ex- 
ponential sensitivity to initial conditions. It is this sensitivity we would like 
to use to distinquish between the cases q = and q > from the previous 
section. 

Consider the so called logistic map which is given by the equation 



e [0, l] • 



The properties of the map depend on the parameter a. If we take, for example, 
a = 3.71, then the Lyapunov exponent is positive, the trajectory is very 
sensitive to the initial value and one has the chaotic behaviour fll2||. It is 



important to notice that if the initial value Xq = 0, then x n = for all n. 




Fig. 1 . Change of x n w.r.t. time n 



It is known [P3] that any classical algorithm can be implemented on quan- 
tum computer. Our stochastic quantum computer will be consisting from 
two blocks. One block is the ordinary quantum computer performing com- 
putations with the output = a/1 — q 2 |0) + q\l). The second block is 
a quantum computer performing computations of the classical logistic map. 
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This two blocks should be connected in such a way that the state first 
be transformed into the density matrix of the form 

p = q 2 P 1 + (1 - q 2 ) P 

where P\ and Po are projectors to the states |1) and |0) . This connection 
is in fact nontrivial and actually it should be considered as the third block. 
One has to notice that P\ and Pq generate an Abelian algebra which can be 
considered as a classical system. In the second block the density matrix p 
above is interpreted as the initial data po f° r the logistic map 

p n+ l = ap n {\ - p n ) 

After one step, the state pi becomes 

p x = aq 2 (l - q 2 )I, 

where / is the identity matrix on C 2 . In paricular, if one has q = then 
Po = -Po an d we obtain p n = P for all n. Otherwise the stochastic dynamics 
leads to the amplification of the small magnitude q in such a way that it can 
be detected. As is seen in Fig.l, we can easily amplify the small q in several 
steps, i.e., within about ten times measurements as in Shor's algorithm. The 
transition from p n to p n+ \ is nonlinear and can be considered as a discrete 
Heisenberg evolution of the variable x n . 

One can think about various possible implementations of the idea of using 
chaotic dynamics for computations. Below we discuss how one can realize 
nonlinear quantum gates on atomic quantum computer. 

5 Atomic Quantum Computer 

Many current proposals for the realization of quantum computer such as 
NMR, quantum dots and trapped ions are based on the using of an atom or an 
ion as one qubit, see |T7|, [TS|, |TJ, |20| . In these proposals a quantum computer 
consists from several atoms, and the coupling between them provides the 
coupling between qubits necessary for a quantum gate. It was proposed 
in [jllj that a single atom can be used as a quantum computer. One can 
implement a single qubit in atom as a one-particle electron state in the self- 
consistent field approximation and multi-qubit states as the corresponding 
multi-electron states represented by the Slater determinant. So, to represent 
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10 qubits one can use an atom with 10 electrons and to represent 50 qubits 
one has to control only around 50 levels in an atom with 50 electrons. 

A possible realization of the standard quantum gates in the atomic quan- 



tum computer by using the electron spin resonance has been discussed in [11 
In this paper we propose that in the atomic quantum computer one can build 
also nonlinear quantum gates because the dynamics of the multi-electron 
atom in the very good approximation is described by nonlinear Hartree-Fock 



equations. Therefore it follows from [10] and the considerations in this paper 



that the atomic quantum computer can solve the SAT problem in polynomial 
time. 

It is well known that in atomic physics the concept of the individual state 
of an electron in an atom is accepted and one proceeds from the Hartree- 



Fock self-consistent field approximation, see for example [pij| . The state 
of an atom is determined by the set of the states of the electrons. Each 
state of the electron is characterized by a definite value of its orbital angular 
momentum I, by the principal quantum number n and by the values of the 
projections of the orbital angular momentum mi and of the spin m s on the z- 
axis. In the Hartree-Fock central field approximation the energy of an atom 
is completely determined by the assignment of the electron configuration, 
i.e., by the assignment of the values of n and / for all the electrons. 

The tensor product structure of states is very important for computations. 
Fortunately a multielectron atom admits such a structure. More exactly, 
instead of the standard tensor product used in quantum computing we have 
to use the Slater determinant to take into account the Fermi statistics. The 
standard computational basis in quantum computing does not have Bose 
or Fermi symmetry. In the atomic case we have to make an appropriate 
modification of quantum gates to take into account Fermi statistics and this 
leads to a new type of entanglement related with Fermi statistics. 

An application of the electron spin resonance (ESR) to process the in- 
formation encoded in the hyperfine splitting of atomic energy levels and to 
build standard linear quantum gates has been considered in [p]] . In this pa- 
per we suggest that in atomic quantum computer one can build also nonlinear 
quantum gates described by the Hartree-Fock equations. 

The Hamiltonian for the N— particle system has the form 

i=l ^ 1 ' i<j 
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In the Hartree-Fock method one takes the iV— particle wave function in the 
form of the Slater determinant 

#(t,n,...,rjv) = Antisym(® 1 (t,r l )...3> N (t,r N )) 

Here the one-particle wave functions rj) satisfy the nonlinear Hartree- 
Fock equations which have the form of nonlinear Schrodinger equation 

i^!4p2 = *(•).,(«, 

where 

r) = f-^- + v (r)+U i (t,r)] $,(t,r)- / dr'W(i, r', r)$;(t, r') 



and 



W (t, r', r) = ^ $* (t, r') V (r', r) ^ (t, r) 



If we consider only the spin dependent part of the wave function of the 
one-electron state 

<Pi (t) 



<P = 

then one can get the nonlinear equation of the form 



Here A is a 2 x 2 matrix and the matrix B depends on (p. By using this 
equation one can describe nonlinear quantum gate. The nonlinearity can be 
tuned by means of magnetic field. 



6 Conclusion 

The complexity of the quantum algprithm for the SAT problem has been 



considered in [10 where it was shown that one can build the unitary matrix 
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Uf in the polynomial time. We have also to consider the number of steps in 
the classical algorithm for the logistic map performed on quantum computer. 
It is the probabilistic part of the construction and one has to repeat compu- 
tations several times to be able to distingish the cases q = and q > 0. Thus 
it seems that the chaotic quantum computer can solve the SAT problem in 
polynominal time. 

In conclusion, in this paper the chaotic quantum computer is proposed. 
It combines the ordinary quantum computer with quantum chaotic dynamics 
amplifier which can be implemented by using the atomic quantum computer. 
We argued that such a device can be powerful enough to solve the NP- 
complete problem in polynomial time. 



References 

M. Garey and D. Johnson, Computers and Intractability - a guide to the 
theory of NP- completeness, Freeman, 1979. 

PW. Shor, Algorithm for quantum computation : Discrete logarithm and 
factoring algorithm, Proceedings of the 35th Annual IEEE Symposium 
on Foundation of Computer Science, pp. 124-134, 1994. 

E. Bernstein and U. Vazirani, Quantum Complexity Theory, in: Proc. 
of the 25th Annual ACM Symposium on Theory of Comuting, (ACM 
Press, New York,1993), pp.11-20. 

C. H. Bennett, E. Bernstein, G. Brassard, U. Vazirani, 
Strengths and Weaknesses of Quantum Computing, |quant-ph / 970 1 0(JT 



L. Fortnow and J. Rogers, Complexity Limitations on Quantum Com- 
putation, ^s.CC/98110231 . 



R. Cleve, An Introduction to Quantum Complexity Theory, |quant- 
ph/9906lll| . 



E. Hemaspaandra, L.A. Hemaspaandra and M. Zimand, 

Almost- Everywhere Superiority for Quantum Polynomial Time, |quant- 

ph/9910033| . 



M. Ohya, Mathematical Foundation of Quantum Computer, Maruzen 
Publ. Company, 1998 



10 



[9] D. S. Abrams and S. Lloyd, Nonlinear quantum mechanics implies 
polynomial-time solution for NP-complete and #P problems, |quant-| 
ph/980104Tl . 



[10] M. Ohya and N. Masuda, NP problem in Quantum Algorithm, |quant 



ph/9809075 



[11] I.V. Volovich, Atomic Quantum Computer, |quant-ph/9911062 . 



[12] M. Ohya, Complexities and Their Applications to Characterization of 
Chaos, Int. Journ. of Theoret. Physics, 37 (1998) 495. 

[13] S.A. Gardiner, J.I. Cirac and P. Zoller, Phys. Rev. Lett. 79(1997) 4790. 

[14] R. Schack, Phys. Rev. A57 (1998) 1634; T. Brun and R. Schack, |quant- 



ph/9807050| . 

[15] I. Kim and G. Mahler, Quantum Chaos in Quantum Turing Machine, 
|quant-ph/9910068| . 

[16] D. Deutsch, Quantum theory, the Church- Turing principle and the uni- 
versal quantum computer, Proc. of Royal Society of London series A, 
400, pp.97-117, 1985. 

[17] J.I. Cirac and P. Zoller, Phys. Rev. Lett., 74 (1995) 74. 

[18] N.A. Gershenfeld and I.L. Chuang, Science, 275 (1997) 350. 

[19] G. Burkard, D. Loss and D.P DiVincenzo, |cond-mat/9808026| . 



[20] A. Ekert and R. Jozsa, Quantum computation and Shor's factoring al- 
gorithm, Reviews of Modern Physics, 68 No. 3,pp. 733-753, 1996. 

[21] I.I. Sobelman, Atomic Spectra and Radiative Transitions, Springer- 
Verlag, 1991. 

[22] Accardi, L. and Ohya, M.: Teleportation of general quantum states, 
|quant-ph/9912087| . 

[23] Fichtner, K.-H. and Ohya, M.: Quantum Teleportation with Entangled 
States given by Beam Splittings, |quant-ph /99 1 2083| . 



11 



